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Summary — A cancellation technique permits measurement of the 
frequency of a source moving relative to an observer without the ob- 
scuring effect of first-order Doppler shifts. The application of this 
method to a gravitational red shift experiment involving the use of 
an earth satellite containing a highly stable oscillator is described. 
The rapidity with which a measurement can be made permits the 
taking of data at various altitudes in a given elliptical orbit. Tropo- 
spheric and ionospheric effects upon the accuracy of results are esti- 
mated. 

Introduction 

ACCORDING to the General Theory of Relativity, 
the time kept by a “proper” clock (or the fre- 
quency of a standard oscillator) is affected by the 
motion of the clock and by its local gravitational poten- 
tial. The latter effect, known as the “gravitational red 
shift,” is expected to result in variations of only 7 parts 
in 10 10 in the rates of clocks in the vicinity of the earth. 
Attempts to obtain experimental verification have been 
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stimulated by the development of highly stable fre- 
quency standards and the availability of earth satellites 
in orbits that provide a significant variation in gravita- 
tional potential. 

In 1956, Singer 1 proposed the measurement of gravi- 
tational red shift by an experiment involving compari- 
son of the time kept by a stable clock in a satellite with 
that of an identical clock on earth. Since then, several 
versions of his experiment have been suggested. Basi- 
cally, all of them involve the use of an atomic clock, or 
equivalent standard, from which, by frequency division, 
a series of precise timing pulses can be derived. These 
pulses would be transmitted as a modulated carrier over 
a one-way link from a satellite to earth, where they 
would be compared with a similar, locally generated 
timing waveform to determine the frequency difference. 
The frequency of the received pulses is directly in- 
fluenced by the rate-of-change of path lengths (first- 
order or radial Doppler shift). This effect is typically 
four or five orders of magnitude larger than the gravita- 

1 S. F. Singer, “Application of an artificial satellite to the measure- 
ment of the general relativistic ‘red shift’,” Phys. Rev., vol. 104, pp. 
11-14; October, 1956. 
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tional red shift. Consequently, data must be^ acquired 
over a sufficiently long time interval (probably several 
days) to insure that path-length measurement errors 
are negligible. At the completion of this experiment, a 
single value is obtained for the average frequency shift 
over many orbital periods. Thus it is not possible by 
this technique to determine the effect of altitude 
variations. 

In order to obtain the altitude dependence of gravita- 
tional red shift, the measurement interval must be 
shortened to a small fraction of the orbital period. How- 
ever, path-length variations make the performance of a 
short-time measurement in a one-way propagation ex- 
periment virtually impossible. But two-way propaga- 
tion experiments offer definite possibilities. Engineers 
are quite familiar with the doubling of the first-order 
Doppler shift, which occurs in transmission of signals to 
and from a radar target. If it were possible to reverse the 
sense ( i.e ., direction) of one of these frequency shifts by 
a mixing operation, then essentially complete cancella- 
tion of first-order Doppler shift could be effected. In 
the short-time measurement system 2 that will be de- 
scribed, such a technique permits data for a single alti- 
tude to be obtained in an interval of 60 seconds or less. 
The following significant advantages result. 

1) By employing an elliptical orbit, frequency dif- 
ferences at various altitudes can be observed, and 
hence the variation of frequency with gravitational 
potential can be obtained. 

2) If the measurement is repeated when the satellite 
is at a particular altitude ( e.g at perigee), a slow 
frequency drift between satellite and earth clocks 
can be determined and removed. Since the quan- 
tity of interest is the variation in frequency dif- 
ference as a function of altitude, a fixed offset fre- 
quency between clocks is inconsequential. There- 
fore, quartz-crystal oscillators can probably be 
employed. 

3) Continuous integration or counting of pulses over 
many orbital periods is not required. Skipped 
cycles or temporary failures would not invalidate 
the experimental results. 

The Short-Time Measurement Principle 



Fig. 1 — Geometry of the measurement problem. 

1) A transmitter (located at point 0) of frequency / 
transmits to the satellite. 

2) The frequency of the signal received at the satellite 
is 


1 — 0 cos a 

/'=/ > ( 1 ) 

(1 -/ 5 2 + 2<#>) 1 ' 2 

where /3 is the magnitude of the satellite velocity 
(measured in the coordinate system of the observer 
at 0) at the instant of reception, normalized in 
terms of the velocity of light; a is the angle (meas- 
ured in the same coordinate system) between the 
propagation path O-S and the velocity vector - /? ; 
<t> is the (positive) difference in gravitational po- 
tential between points 5 and 0, normalized in 
terms of the square of the velocity of light. 

3) The satellite contains an oscillator of frequency 2/ 
(measured on the satellite). The received signal 
frequency /' is subtracted from 2/ by mixing 
within the satellite to yield a frequency difference 

f - 2 !/-/'. ( 2 ) 

4) The mixture signal/" is transmitted back to earth. 
(Reception and retransmission at the satellite are 
assumed to be virtually instantaneous.) 

The frequency of the received signal is measured as 


r =/' 


(i - £ 2 + 2 4>y n 


1 + /3 cos a 
Substitution for f" in (3) yields 


T 1-/3 cos a ~ir(l -/3 2 + 20) ,/ H 

/"' = / 2 . 

L (1 - /3 2 + 2<£) 1/2 J L 1 + /8 cos a J 


(3) 


(4) 


The principle employed in a short-time measurement 
will be illustrated by a simplified analysis. In the interest 
of simplicity, the measurements are assumed to take 
place at an earth pole (0 in Fig. 1). (An outline of the 
derivation of the following expressions is given in the 
Appendix.) 

2 R. S. Badessa, R. L. Kent, and J. C. Nowell, “Short-time meas- 
urement of time dilation in an earth satellite,” Phys. Rev . Lett., 
vol. 3, pp. 79-80; July, 1959. 


By means of a power-series expansion, retaining second- 
order terms in /3 and first-order terms in </>, we obtain 

/"' « /(I + 2* - 0 2 ). 

The first-order Doppler term has been cancelled, and the 
received signal frequency differs from the ground trans- 
mitter frequency by 

A / « f(2<t> - 0 2 ). 


(5) 
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The two terms in (5) are apparently inseparable experi- 
mentally. The first term, 2 $/, is the desired gravitational 
frequency shift. The second term, — /3 2 /, is the well- 
known second-order, or transverse, Doppler shift of 
special relativity. This term, which has been verified 
experimentally, can be calculated easily and subtracted 
from the measured frequency difference to obtain the 
residual gravitational effect. 

The frequency difference in the general case of an ob- 
server not located at the pole (Fig. 2) is given by 

A/ = f[2<{> - (ft 2 — /3 0 2 ) - 2ftft cos (a } + 72 ) 

— ft> 2 (l + r sin 72 )], ( 6 a) 

where fi 3 , /3 0 , a 2 , and 72 are as shown in Fig. 2, measured 
in a nonrotating coordinate system with origin at the 
earth’s center and r is the range, normalized in terms of 
the earth’s radius. The term in ( 6 a) containing the prod- 
uct ftol3 9 is of the order of 10 -10 , and cannot be neglected 
in the present experiment. This term must be retained, 
and appropriate values substituted. Present orbital 
computational accuracy is more than adequate for this 
purpose. 

Satellite Equipment 

As we have indicated, the operations to be performed 
in the satellite are as follows: 

1 ) Amplify a received signal of frequency/'. 

2) Mix the received signal with that of the satellite 
“clock” to obtain 2 /—/'. 

3) Amplify and transmit a signal of frequency 2 /—/'. 

In order for the cancellation of first-order Doppler 
shifts to be effective, the signals received and trans- 
mitted at the satellite must be very nearly of the same 
frequency. One way of avoiding the ensuing “spillover” 
problem is to employ pulsed waveforms. The received 
signal pulses, 1-10 jusec long, are stored briefly in a delay 
line or equivalent device. After the reception of each 
pulse, the receiver is turned off and the transmitter 
turned on. Fig. 3 is a block diagram of the satellite 
equipment required to instrument this technique. This 
configuration is intended only for purposes of illustra- 
tion. An operating frequency range of 500-1000 me, 
which is expected to fall between regions of severe iono- 
spheric and tropospheric propagation anomalies, has 
been tentatively selected. A factor in the final choice of 
an operating frequency is the desirability of employing 
transistor circuitry. For ranges up to 10,000 miles, an 
RF output of approximately 20 mw (average), 0.4 watt 
(peak) is required from the satellite transmitter. 

Although duplication of units is not indicated in the 
block diagram, it is particularly desirable to include 
more than one crystal frequency standard in the satel- 
lite, to permit independent measurements and to im- 
prove system reliability. The use of three crystal oscil- 
lators in a common temperature-stabilized oven is sug- 
gested. The crystals, which could be offset from one 
another in frequency for identification purposes, might 



Fig. 2 — Geometry of the two-way propagation experiment. 



Fig. 3 — Illustrative satellite-equipment configuration. 


be selected with different known temperature coeffi- 
cients. Then variations in satellite temperature, which 
could possess a deceiving similarity to the orbital period 
of the satellite itself, will affect the three oscillators in 
varying degrees. Thus, a means for separating relativ- 
istic frequency shifts from temperature-dependent ef- 
fects is provided. Preliminary power and weight esti- 
mates for the satellite equipment, including the weight 
estimate for batteries, are 3 watts and 35 pounds. 

The stability of the crystal oscillators is of paramount 
importance and will be discussed separately. The over- 
all phase shift through the receiver and transmitter 
should remain stable within 30° over the measurement 
interval 20-60 seconds, in order not to degrade the 
measurement accuracy which should be determined es- 
sentially by the short-term stability of the crystal os- 
cillators. 

Ground Equipment 

An illustrative block diagram of the major portion of 
the required ground equipment is shown in Fig. 4. It 
functions as follows: 

1 ) A signal of frequency / is transmitted. 

2) A signal of frequency /'" is received from the satel- 
lite and compared in frequency with / by means of 
a mixer and a synchronous detector. 

3) The frequency-difference signal is recorded on tape 
together with a fixed-frequency reference signal. 
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The transmitted and received signals both consist of 
series of pulses (each 1-10 ^sec long). Time-delay dis- 
crimination is provided in the receiver for the rejection 
of transmitter feed-through signals and other undesired 
signals. In addition, a frequency offset of 1 part in 10 6 
is purposely introduced between the ground and satel- 
lite frequency standards. This insures a small frequency 
separation between the signals transmitted and received 
on the ground, thereby facilitating resolution of the 
satellite signal from clutter by means of narrow-band 
filtering, without materially affecting first-order Doppler 
cancellation. 

The coherent waveform and narrow-band filtering 
( ± 1 cps) employed in the ground receiver make possible 
the low requirement on satellite transmitter power (20 
mw, average). Doppler tracking is unnecessary because 
the frequency of the signal received on the ground is 
virtually independent of first-order Doppler shift. How- 
ever, the first-order Doppler shift is present in the signal 
received on the satellite. For simplicity, broad-band 
amplification, rather than Doppler tracking, is em- 
ployed there, at the expense of increased ground- 
transmitter power. For a 10 meter 2 antenna on the 
ground, the power required to yield a clean signal on the 
satellite is approximately 500 watts (average), 10 kw 
(peak). 

The ground station requires means for storage, as 
well as immediate processing, of the receiver output 
signal. The desired information is contained in the small 
variations in output frequency with orbital position. 
Since orbital position will be known as a function of 
time, a recording of the output signal together with a 
time signal will suffice. 

Processing of the data involves measurement of a 
small frequency difference between two signals. This is 
best accomplished by measuring the total phase dif- 
ference accumulated in a specific time interval. For ex- 
ample, a frequency difference of 0.1 cps in a 60-second 
interval results in a total accumulated phase difference 
of 2160°. This can be measured precisely by manual 
phase-tracking techniques, with the use of an oscillo- 
scope as a visual tracking indicator. 

Oscillator Stability 

The frequency standards used in this experiment must 
meet two separate stability requirements. First, all 
random and periodic frequency excursions must be at 
least an order of magnitude smaller than the total effect 
that is to be measured. Second, any long-term frequency 
drift in one orbital period should not exceed 20 to 30 per 
cent of the total desired effect. Perhaps the reasons for 
these stability requirements can best be illustrated by 
plotting a set of hypothetical data that might be ob- 
tained in this experiment. Fig. 5 shows a plot of the 
hypothetical frequency difference against time. For 
simplicity, it is assumed that data are acquired alter- 
nately at apogee and at perigee for several successive 
sightings. The clusters of points at times A, B, C, and so 



Fig. 4 — Illustrative ground-equipment configuration. 



Fig. 5 — Hypothetical data to illustrate oscillator 
stability requirements. 

on, represent several readings of frequency difference ob- 
tained during the intervals when the satellite is visible. 
The frequency difference between the two dashed lines 
for apogee and perigee readings is expected to be ap- 
proximately 4 parts in 10 10 . Therefore, the undesired 
random and periodic frequency variations of the oscil- 
lators that are indicated by the dispersion of the clusters 
must be no greater than a few parts in 10 n . 

The slope of either of the dashed lines in Fig. 5 repre- 
sents the difference between the long-term drifts of the 
oscillators, and it is clear that a monotonic drift per or- 
bital period of from 20 to 30 per cent of the total effect 
can easily be tolerated. For the proposed experiment, 
a drift of 1 part in 10 10 in a few hours, or approximately 
1 part in 10 9 per day, is permissible. A small fixed offset 
frequency between the oscillators would merely shift the 
base line in Fig. 5, and therefore would have no effect 
upon the desired measurement. 

Actual measurements taken on commercially avail- 
able quartz-crystal oscillators indicate that these oscil- 
lators essentially meet the stability requirements. Table 
I shows some data obtained in a comparison test. The 
1-mc output signals of two crystal oscillators were multi- 
plied up and mixed down in such a way as to magnify 
their frequency difference by a factor of 1000. The re- 
sulting beat frequency was then determined by means 
of a counter operating over periods of 10 and 100 
seconds. The last digit of the counter output represents 
parts in 10 l ° for the 10-second count, and parts in 10 u 
for the 100-second count. In either case, the counter has 
a starting error of one count. Thus, the differences in- 
dicated between the two oscillators may be entirely ac- 
counted for by the counter error. The oscillator error is 
shown to be no greater than 1 part in 10 11 when it is 
averaged over successive 100-second intervals. 










762 


PROCEEDINGS OF THE IRE 


TABLE I 

Frequency Comparison Data for Two Hermes Model 101c 
1-mc Crystal Frequency Standards 


10-Second Counting Interval; 

Total Elapsed Time, 
Approximately 45 Minutes. 

100-Second Counting Interval; 
Total Elapsed Time, 
Approximately 15 Minutes 

93101 

31005 

93100 

31006 

93100 

31006 

93100 

31005 

93100 

31006 

. 

31006 

; 

31006 

93100 


93100 


93101 


93100 


93101 



The oscillators that were tested employed vacuum- 
tube circuits; data are not available for a transistor 
version that is now being developed. However, the sta- 
bility is believed to be determined primarily by the 
crystal and its oven, and not by the oscillator circuitry. 
Therefore no degradation in stability is expected to re- 
sult from a change to transistors. 

In order to transfer the oscillator stability to a trans- 
mitter operating in the UHF range, a chain of frequency 
multipliers is required. Two methods of frequency multi- 
plication, both of which use variable-capacitance diodes 
as the multiplying elements, are being considered. 
Leeson and Weinreb, 3 in a detailed study of low-order 
frequency multipliers with nonlinear reactances, demon- 
strated that frequency doublers with only 1-db loss 
could be constructed. A chain of doublers that are now 
being developed to go from the oscillator frequency to 
the transmitter frequency is expected to yield a loss of 
only 10-15 db. The second approach to the frequency- 
multiplication problems is the use of one or two large 
frequency jumps. However, experimental results thus 
far indicate that in multiplying by a factor of 100, the 
loss is approximately 40 db. 

Propagation Effects 

As we have stated, the success of the described meas- 
urement depends upon the cancellation of the first-order 
Doppler shift, an effect that is four or five orders of 
magnitude greater than the gravitational red shift. For 
adequate cancellation, it is essential that random fluc- 
tuations in the effective path lengths upward and down- 
ward be small, or that they be well correlated. Phase 
fluctuations along the propagation paths are known to 
be caused by turbulent irregularities in the troposphere 
and in the F 2 region of the ionosphere. Extensive dis- 

8 D. B. Leeson and S. Weinreb, “Frequency multiplication with 

nonlinear capacitors — a circuit analysis,” Proc. IRE, vol. 47, p 

2076-2084; December, 1959. 


April 


cussions of these and related effects have been pub- 
lished. 4 " 8 * 10 

It is desirable to select an operating frequency in a 
region in which the total rms phase error that is ex- 
pected to arise from propagation anomalies falls well 
within the total permissible phase error of the experi- 
ment. If the frequency shift is computed by observing 
the total accumulated phase difference in a specified 
time interval T , then the resulting frequency error can- 
not exceed A\p/i rT cps, where A\p is the total uncancelled 
phase fluctuation. At 500 me, an rms phase error of 0.5 
radian results in a 1 per cent error in frequency for a 
60-second measurement interval. Such an error would 
be acceptable, since it is less than the value of the ex- 
pected oscillator instability. 

Tropospheric fine structure is measured by variations 
in the dielectric constant. The rms one-way phase 
fluctuation obtained by Muchmore and Wheelon 4 is 
given by 

(2/ 0 Z) 1/2 7rAe 

Tprma — > for L » / 0 , (6b) 

X 

where /o represents the scale of the irregularities (prob- 
ably 200 feet) ; L is the propagation path length through 
the turbulent medium (taken to be 20,000 feet) ; Ae is 
the rms deviation of the dielectric constant from unity 
(probably 10" 6 ) ; and X is the wavelength. At a frequency 
of 1000 me, the one-way phase fluctuation is computed 
as approximately 0.6° rms. This value is believed to be 
representative for a turbulent mass of (optically) clear 
air. In the presence of rain clouds, the interiors of which 
are very turbulent and contain large fluctuations in di- 
electric constant, the situation is somewhat different. 
If we take as representative values for a rain cloud 
lo = 20 feet, Z, = 5000 feet, Ae = 20X10~ 6 , the one-way 
fluctuation in phase at 1000 me is found to be approxi- 
mately 1.6° rms. Thus, a net phase fluctuation of 2° or 
3° rms might be typical for 1000-mc one-way propaga- 
tion through the troposphere on a cloudy day. 

The tropospheric phase-fluctuation measurements of 
Herbstreit and Thompson, 10 at 1046 me, appear to be in 

4 R. B. Muchmore and A. D. Wheelon, “Line-of-sight propagation 
phenomena. I. Ray treatment,” Proc. IRE, vol. 43, pp. 1437-1449; 
October, 1955. 

6 J. L. Pawsey and R. N. Bracewell, “Radio Astronomy,” Claren- 
don Press, Oxford, Eng., pp. 347-351; 1955. 

8 H. G. Booker, “Tne use of radio stars to study irregular refrac- 
tion of radio waves in the ionosphere,” Proc. IRE, vol. 46, pp. 298- 
314; January, 1958. 

7 D. F. Martyn, “The normal F region of the ionosphere,” Proc. 
IRE, vol. 47, pp. 147-155; February, 1959. 

8 C. O. Hines, “Motion in the ionosphere,” Proc. IRE, vol. 47, 
pp. 176-186; February, 1959. 

9 “Propagation Factors Affecting Long-Range UHF Radars at 
High Latitudes,” Defence Res. Telecommun. Establ., Ottawa, Ont., 
Can., Rept. No. 41-1-3; December, 1958. 

10 J. W. Herbstreit and M. E. Thompson, “Measurements of the 
phase of radio waves received over transmission paths with electrical 
lengths varying as a result of atmospheric turbulence,” Proc. IRE, 
vol. 43, pp. 1391-1401; October, 1955. 
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general agreement with theory. The standard deviation 
of the phase shift over a 10-mile tropospheric path ex- 
ceeded 3° fifty per cent of the time, and 5.5° ten per cent 
of the time. Moreover, the fluctuation period was found 
to be approximately 1 hour for large fluctuations, with 
minute-to-minute variations in phase shift of only 0.3° 
rms. Since the tropospheric phase fluctuation, from (6b), 
varies directly with the operating frequency, we con- 
clude that the contribution of the troposphere to meas- 
urement errors in this experiment is negligible. 

Ionospheric fluctuations are measured by variations 
in the free-electron density. Turbulent irregularities in 
the F 2 region are maximum near midnight (local time) 
and occur in a region roughly 55 km thick at a nighttime 
altitude of approximately 300 km. The one-way rms 
phase fluctuation is shown by Booker 6 to be 

//*otf \ 1/2 

* rms = i — ) r e \AN , (7) 

where h 0 , H are analogous to / 0 , L for the troposphere 
(taken here to be 1 km and 55 km, respectively); r c is 
the classical radius of the electron, 2.8 X10 -16 meter; 
and A N is the rms fluctuation in free-electron density. 
A value of A N~ 10 u /meter 3 is believed to be representa- 
tive of the severe fluctuations encountered during peri- 
ods of intense sunspot activity. Use of this value yields 
a phase fluctuation of 0.44 radian rms at 1000 me. The 
ionospheric effect, from (7), varies inversely with the 
operating frequency. The desired relativistic effects, 
from (5), vary directly with frequency. Therefore, the 
fractional measurement error caused by ionospheric 
turbulence will vary inversely as the square of the oper- 
ating frequency. Recent UHF measurements indicate 
that the phase fluctuation in propagation through the 
ionosphere is much less than one radian rms during a 
daytime period of normal ionospheric activity. 

The Doppler-cancelling property of the two-way 
propagation experiment will effect cancellation of phase 
(or frequency) fluctuations resulting from propagation 
anomalies, as well as the radial Doppler shift resulting 
from the motion of the satellite, to the extent that cor- 
relation exists between the fluctuations introduced in 
the upward and downward traversals of the turbulent 
medium. The correlation coefficient is still to be esti- 
mated, in order to determine the residue of the calcu- 
lated one-way phase fluctuation that will remain un- 
cancelled in the experiment. It is known that the iono- 
sphere rotates with the earth, and hence the upward- 
and downward-ray paths pass through essentially the 
same region of the ionosphere. The round-trip time de- 
lay, r, associated with the largest anticipated range 
(10,000 miles) is 0.1 second. The value of the correlation 
coefficient for phase fluctuations displaced up to 0.1 
second in time is therefore required. Experimental 
data 6 9 indicate maximum fluctuation rates of approxi- 


mately 4 cycles per minute, occurring during periods of 
magnetic activity. Assuming an exponential correlation 
function, we have 

p( T ) = (g) 

where co 0 = 27r/15 =0.42 radian/second. 

The residual (uncancelled) phase error is then given 
by 

tyrms = *rn.s[2(l - p)]«* (9) 

which, for small values of co 0 r, is approximately 

A^rms ~ ^rm 8 (2cd0T) 1/2 , (o) 0 T « 1). (10) 

For r = 0.1 second, the quantity (24O 0 r) 1/2 in (10) is equal 
to 1/3.5. The Doppler-cancellation process is therefore 
effective in reducing ionospheric phase fluctuations by 
a factor of 3.5, or better. The uncancelled phase shift 
is found to be 7.4° rms at 1000 me, or 15° rms at 500 me, 
for a path length (one-way) of 10,000 miles. Even at the 
lower operating frequency, the expected error for a 60- 
second measurement interval is only 0.0014 cps, or 0.4 
per cent of the expected maximum red shift. 

We conclude that for operating frequencies in the 
range 500-1000 me, propagation anomalies will produce 
a negligible deterioration in the accuracy of the 
measurement. 

Appendix 

We shall outline a method for computation of the 
Doppler frequency shift associated with the propagation 
of electromagnetic signals between points in nonuniform 
motion in a gravitational field. An approximate result is 
derived here; an exact derivation has been made by 
Friedman. 11 

If a gravitational field in space is described in a sys- 
tem of coordinates x u oc 2 , x z , x 4 , where x± — ct, and / 
represents “coordinate time,” then the relation 

ds 2 = H gurdx^dx, (11) 

M=1 v=l 

holds between each locally measured interval ds and the 
corresponding coordinate differences, dx For events 
occurring at a single point, the locally measured interval 
is given by 

ds = cdr> (12) 

where dr denotes “proper time” at that point. The 
quantities g „„ define the metric properties of space ( i.e ., 
the gravitational field). For a spherically symmetric dis- 
tribution of mass about the origin (Fig. 6), the g are 
given by the Schwarzschild solution to the field equa- 
tions of general relativity: 

11 M. H. Friedman, Dept, of Physics, Cambridge, Mass., 

private communication; October, 1959. 
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Fig. 6 — Geometry of propagation between two points moving in 
a spherically symmetric gravitational field. 


£44 = 1 + 20 


— gftA ^ 

+ 


20 

1 + 20 


A jj.Yj, 


M = 1,2,3 
^ — 1,2, 3, 


(13) 


where =1 (ju = y) ; 5 M „ = 0(^u v) ; r is the distance to the 
origin from the point x u *2, .r 3 ; and 0 is the normalized 
gravitational potential at that point. 

Consider a transmitter (Fig. 6) at point 1 in gravita- 
tional potential 0, whose frequency / is precisely con- 
trolled by a resonance of atoms, traveling with respect 
to the origin with (normalized) velocity ft. The fre- 
quency transmitted may be regarded as some number 
of events per unit time. As viewed in the jc m coordinate 
system, this frequency is given by 


dri 



dn dr 1 dr\ 

dr 1 dt\ dli 


(14) 


Substitution for g^ v and ds in (11) yields 

— = [1 + 24, - j8i*(l - «icos**0l 1 '*, (15) 

dh 


in which €1 = (20i)/(l +200 ; ft is the angle between ft 
and the radius to point 1. The observer at point 2 is as- 
sumed to utilize a frequency standard of like constitu- 


tion to that at point 1. By the same method, the fre- 
quency received at point 2 (viewed in the coordinate 
system) can be related to the frequency measured. Thus 

dr 2 

/*=/'— ( 16 ) 
dt 2 


and 


— = [1 + 202 - 02 2 (1 - €2 COS” 6 .) ]*/*, (1 7) 

dt 2 


where dt 2 is the (coordinate) time of arrival of the ray at 
point 2. Since the number of events is preserved, 

dr 2 dr\ 

dn = /' dt 2 = f dt\. (18) 

di 2 dh 

Therefore 


/ =/ 


dh 

dh 


1 + 2<f>, - /3, 2 (1 - €1 cos 2 0i)"| I/2 

1 + 202 — ft 2 (l — €2 COS 2 $ 2 ) - 


(19) 


The relationship between the coordinate times of recep- 
tion and transmission is intractable because the ray 
velocity varies slightly along the path, and the path it- 
self is curved. However, it can be shown that neglect of 
these effects leads to an error of the order of the satellite 
velocity multiplied by the gravitational potential, 
which is only approximately 10~ 14 for earth satellites. 
For present purposes, the derivation can be completed 
by a classical computation of the rate-of-change of path 
length: 

dh 1 — ft cos a 2 

« ( 20 ) 

dU 1 + ft cos a.\ 


Moreover, ti and e 2 in (19) can be neglected. Thus, 
^ T I + 20x - ft 2 T' 2 / l -ft cos aA 

f / Ll + 20 2 -ft 2 J \l+ft C osai / 


( 21 ) 


Acknowledgment 

The authors wish to acknowledge the assistance given 
by Prof. M. H. Friedman, Prof. J. W. Graham, and V. J. 
Bates. They would also like to express their apprecia- 
tion for the invaluable guidance of Prof. J. R. Zacharias. 



